PLOS Genetics

Meiotic recombination is a key driver of evolution in sexually reproducing organisms, reshaping genetic diversity by generating novel allelic combinations. The rate of recombination varies substantially across living organisms depending on cis- or trans-acting genetic elements, as seen in many species, including the yeast Saccharomyces cerevisiae. Here, we report on an experimental evolution-based study to better understand the factors shaping this natural variation. Starting with a genetically diverse population of S. cerevisiae, we have carried out recurrent divergent selection on recombination rate using a fluorescence-based sorting approach in four independent lineages. After ten generations, we observed an average response of recombination rate of +28% after positive selection and -24% after negative selection, within the interval used for selection. In the adjacent region, however, we observed a weaker response in the opposite direction, and no response in four other unlinked genomic regions. Whole-genome sequencing of individuals selected for high recombination revealed mixed outcomes in the four independently evolved lineages. All four lineages showed selection for high recombination locally, with particular haplotypes heavily favored and sequence- or structural variation-based heterozygosity selected against within the selection interval. However, only two of the four lineages showed increases in genome-wide recombination rate. Overall, this experimental evolution approach provides original and useful insights into the evolvability of the meiotic recombination rate and the associated genetic determinants.

Translation reinitiation (REI) is one of the most important gene-specific regulatory mechanisms by which eukaryotic cells influence expression of main translons, for example during highly conserved integrated stress response (ISR). In S. cerevisiae, expression of the key stress response gene, GCN4, is controlled by an intricate interplay among four short upstream translons (uTranslons, formerly uORFs), resulting in high or low levels of REI at GCN4 depending on the growth conditions. Under nutrient rich conditions, GCN4 expression is repressed, but upon amino acid starvation, it is derepressed, despite of a general translational shut down. Capitalizing on our screening reporter system, we identified three new factors influencing efficiency of REI after translation of GCN4 uTranslons: Rai1p (an RNA quality control and processing factor), and Ssz1p and Zuo1p (members of the Ribosome Associated Complex [RAC]). Importantly, we showed that depletion of these factors deregulated derepression of Gcn4p synthesis under starvation. Furthermore, we found that similar to RAC, Rai1p associates with 40S subunits and actively translating ribosomes. We also explored interactomes of these proteins. Collectively, we present three previously unknown factors that co-regulate stress response to amino acid starvation in the budding yeast by unique mechanisms.

The CenIR regulatory system of Clostridioides difficile comprises a predicted transcriptional repressor, CenI, and a predicted membrane metalloprotease, CenR. The physiological role of CenIR and activating signal(s) are not known. CenIR belongs to the BlaIR family of regulators that mediate resistance to {beta}-lactam antibiotics. In canonical BlaIR systems, binding of a {beta}-lactam to the extracellular transpeptidase domain of BlaR triggers proteolysis of BlaI and thus induction of a closely linked {beta}-lactamase gene. However, CenR lacks a {beta}-lactam-binding domain and transposon mutagenesis indicated CenI is essential for viability even when {beta}-lactams are not present. Here we confirmed essentiality of CenIR and determined its regulon contains [~]12 genes, including an exported protein of unknown function (CDR_0474) that is induced about 500-fold and a peptidoglycan hydrolase (Cwp6) that is induced about 7-fold when cells are depleted of CenIR. There are no essential genes or {beta}-lactamases in the regulon. Phenotypic characterization of CenIR-depletion strains revealed slower growth, mild elongation and cell lysis. Deletion of cdr_0474 corrected all three defects, while deletion of cwp6 only rescued the lysis phenotype. It was possible to delete cenIR in either a {Delta}cdr_0474 or {Delta}cwp6 background. We propose that CenIR is essential because its absence leads to lysis due to Cwp6 overproduction. Bioinformatic analyses revealed the predicted extracellular sensing domains in annotated "BlaR" proteins are diverse. Thus, BlaIR systems are not dedicated to defense against {beta}-lactams but probably enable bacteria to adapt to a variety of environmental stimuli. ImportanceMany of the regulatory systems for controlling cell envelope biogenesis and stress responses have yet to be studied. Here we characterize a Clostridioides difficile BlaIR-like regulatory system that we have named CenIR for cell envelope. Unlike canonical BlaIR systems, which bind {beta}-lactams and induce a {beta}-lactamase, CenIR lacks a {beta}-lactam binding domain and is essential for viability even in the absence of antibiotics. We identified the genes in the regulon and found that CenIR is essential because its absence leads to overproduction of the Cwp6 peptidoglycan hydrolase. We also show that most annotated BlaIR-like systems lack a {beta}-lactam-binding domain, from which we infer that these systems have much broader physiological roles than generally appreciated.

Transcription factor IIIC (TFIIIC) is a multi-subunit protein complex that recruits RNA polymerase III (Pol III) to the majority of its target genes. Evolutionarily conserved overlap between TFIIIC binding sites and the structural maintenance of chromosomes (SMC) complexes suggested a role for TFIIIC in SMC regulation and 3D organization of eukaryotic genomes; but the evidence has remained largely correlational due to the essential role of TFIIIC in RNA Pol III transcription. Here, we directly tested the function of TFIIIC in SMC complex regulation by using auxin inducible depletion in C. elegans. We performed Hi-C and ChIP-seq analyses upon acute depletion of TFTC-3, an essential TFIIIC subunit, and RPC-1, the catalytic subunit of RNA Pol III. Our results show that TFIIIC regulates the localization of two different types of SMC complexes, cohesin and condensin. TFIIIC is also required for increased 3D contacts between distant tRNA genes located on the same or different chromosomes. Depletion of individual SMC complexes did not significantly reduce the intrachromosomal tRNA gene contacts, suggesting redundancy or an independent mechanism mediating these contacts. Together, our study demonstrates an RNA Pol III independent function for TFIIIC, regulating binding of both cohesin and condensin, as well as the 3D organization of tRNA genes.

The calcium-sensing receptor (CaSR) is a class C G protein-coupled receptor (GPCR) with an important role in calcium homeostasis, activating mutations of which cause hypocalcemia. Despite the receptors clinical importance and evidence suggesting {beta}-arrestin-1/2 can modulate CaSR signaling, the residues within CaSR required for {beta}-arrestin recruitment and the role of {beta}-arrestin in desensitization, trafficking and signaling are ill defined. We confirmed that {beta}-arrestin-1 and {beta}-arrestin-2 are recruited to CaSR upon receptor activation. Deletion of the distal cytoplasmic region of CaSR, which replicates several mutations identified in individuals with hypocalcemia, reduced {beta}-arrestin-1/2 recruitment and enhanced receptor signaling. Examination of the receptor cytoplasmic region identified three regions of serine and threonine residues that resemble phosphorylation codes identified in other GPCRs. Mutation of each of these residues to alanine demonstrated one region between amino acids 1003-1011 is important for both {beta}-arrestin-1 and {beta}-arrestin-2 recruitment, while disruption of a second region between residues 976-981 impaired {beta}-arrestin-2 recruitment. Alanine mutagenesis of these residues also reduced CaSR signaling and impaired receptor internalization suggesting an important role in receptor desensitization. Three variants in the ClinVar dataset reported in disorders of calcium homeostasis were identified and two, T1006M and T1008P, were shown to reduce {beta}-arrestin recruitment and enhance CaSR signaling. Further analysis of the T1006M variant showed it reduced CaSR internalization. Thus, we have identified two regions within the CaSR cytoplasmic region that are important for {beta}-arrestin recruitment, receptor desensitization and internalization. Mutation of residues within these regions may represent another mechanism for the pathogenesis of hypocalcemia.

LIN-12/Notch signaling regulates C. elegans vulval development via cell fate specifications in the gonad and epidermis. In the somatic gonad LIN-12/Notch activity specifies the anchor cell (AC) versus ventral uterine cell (VU) fates, with VU receiving more signal. The AC secretes epidermal growth factor (EGF) which induces the underlying vulval precursor cells (VPCs) to adopt vulval fates. In the VPCs the secondary vulval fates are specified by LIN-12/Notch activity. We previously reported that the AGEF-1, an Arf GEF homologous to ArfGEF1 and ArfGEF2, the ARF-1 GTPase, and the adaptor protein complex 1 (AP-1) inhibit LET-23/EGF receptor (EGFR) signaling in the VPCs by antagonizing LET-23/EGFR basolateral localization. Here we report that AGEF-1, ARF-1 and AP-1 regulate LIN-12/Notch signaling during somatic gonad and vulval development. The lin-12(n302) partial gain-of-function causes a potent Vulvaless phenotype due to a lack of AC specification. We demonstrate that loss of AGEF-1, ARF-1 or AP-1 restored the AC fate in lin-12(n302) animals, indicating that AGEF-1/ARF-1/AP-1 promotes LIN-12/Notch signaling in the somatic gonad. Interestingly, loss of AGEF-1, ARF-1 or AP-1 also induced ectopic vulval secondary fates in lin-12(n302) animals, indicating that AGEF-1/ARF-1/AP-1 inhibits LIN-12/Notch in the VPCs. Using a LIN-12/Notch biosensor we demonstrate that loss of UNC-101/AP-1 results in decreased signaling in the VU cell and increased signaling in the VPCs that correspond with decreased expression levels of LIN-12/Notch and LAG-1/DSL ligand in the presumptive AC and VU while also causing increased apical localization of LIN-12/Notch in the VPCs. We hypothesize that the differential regulation of LIN-12/Notch signaling could reflect different trafficking pathways in epithelial cells (VPCs) versus non-epithelial cells (AC and VU). Our results indicate that the AGEF-1/ARF-1/AP-1 trafficking pathway maintains the VPC cell fate patterning by limiting both LET-23/EGFR and LIN-12/Notch signaling. Author summaryCell signaling and membrane trafficking are highly interconnected processes whereby membrane trafficking can regulate signal transduction pathways and vice versa. We previously demonstrated that the ARF-1 GTPase, the downstream AP-1 clathrin adaptor and upstream activator AGEF-1 antagonize the membrane trafficking of the Epidermal Growth Factor Receptor (EGFR) and hence signaling during C. elegans vulva induction. Strong loss of the ARF-1 GTPase pathway resulted in ectopic vulval induction. Here we demonstrate that the ARF-1 GTPase pathway differentially regulates Notch signaling to regulate vulva induction. In the somatic gonad it promotes Notch signaling to regulate the specification of the anchor cell which secretes the inductive signal. In the vulva precursor cells, the ARF-1 GTPase pathway antagonizes Notch signaling which cooperates with EGFR signaling to induce the vulval cell fates. We hypothesize that the differential regulation of Notch signaling by the ARF-1 GTPase pathway could be a result of more complex membrane trafficking pathways in polarized epithelial cells (vulva precursors) versus non-epithelial cells in the developing somatic gonad. Thus, the AGEF-1/ARF-1/AP-1 antagonizes both EGFR and Notch signaling in ensuring that only three of the six vulval precursor cells adopt are induced.

DDI2 and DDI3 (DDI2/3) are duplicated genes in Saccharomyces cerevisiae that exhibit strong induction by a transcription factor Fzf1 in response to chemical treatments like cyanamide (CY) and methyl methanesulfonate (MMS). Although, like DDI2/3, SSU1, YHB1 and YNR064C also contain an Fzf1-binding consensus sequence CS2 and are coordinately regulated by Fzf1, these genes are only modestly induced by CY and MMS. To identify additional cis-acting elements in the DDI2/3 promoter, we made DDI2/3 promoter deletions in a reporter system and identified upstream repressing sequences (URS) spanning 480 nucleotides. To test a hypothesis that the chromatin structure constitutes the URS, we utilized a yeast strain capable of histone H3/H4 depletion by shifting carbon sources. Following histone depletion, DDI2/3 were strongly induced in an Fzf1 dependent manner, while YHB1 was repressed. Interestingly, under histone depletion conditions, CY or MMS treatment further increased expression of all Fzf1-regulated genes to comparable levels in an Fzf1 dependent manner. A genome-wide MNase-seq analysis showed that CY treatment reduced the nucleosome occupancy at the mapped DDI2/3 URS region in wild-type cells, but not in in fzf1{Delta} cells. These findings collectively indicate that Fzf1 plays dual roles in regulating the DDI2/3 response to CY. Firstly, it binds CS2 and serves as a transcription activator. Secondly, it is required for the chromatin remodeling at URS. This two-tier regulation at the DDI2/3 promoter helps to explain why DDI2/3 achieve much higher fold induction by CY and MMS than other Fzf1-regulated genes, suggesting Fzf1 to be a candidate pioneer transcription factor.

Changes in regulatory sequences controlling the timing and activity of gene products underlie much of natural phenotypic variation. Yet, what these changes are and how they impact gene expression remain largely unknown. To address this question, we investigated how transcriptional activity and homeostatic responsiveness of orthologous promoters of the metabolic gene TDH3 evolved among Saccharomyces yeast. We found that promoter expression level increased specifically in the S. cerevisiae lineage and that a substantial part of this increase was caused by genetic variants located between the well-characterized, conserved binding sites for two direct transcriptional regulators. These nucleotide changes altered the promoters expression levels while leaving the expression dynamics conserved. Further, the effects of these nucleotide changes were only seen in the presence of a third transcription factor, TYE7p, which is recruited by the other transcription factors through protein-protein interactions. These results suggest that the cis-regulatory changes act through their influence on the collective assembly/activation of the transcription factors, and that changes acting through such a mechanism can allow distinct parts of gene expression, such as expression level and dynamics, to be tuned separately.

Deciphering the genetic architecture of complex quantitative phenotypes remains challenging in quantitative genetics. These traits not only depend of multiple genetic factors but are also established over time and environments. Although quantitative genetics has investigated the genetic determinism of phenotypic plasticity in contrasted environmental conditions, the time related phenotypic plasticity has received less attention. Here we proposed a multivariate Bayesian framework, the Bayesian Varying Coefficient Model, designed for analysing the genetic architecture of the time related phenotypic plasticity by a multilocus approach. We applied the BVCM to time series phenotypes measured at various time scales (daily, monthly, yearly) across a diverse set of biological species. We included in this study: yeast (Saccharomyces cerevisiae), fungi (Fusarium graminearum), eucalyptus (Eucalyptus urophylla x E. grandis), and sweet cherry tree (Prunus avium). The BVCM results were compared with those obtained with a known genome-wide association method carried out time by time. For all species and traits, the BVCM was able to detect the major QTL identified by marker-trait association methods and revealed additional genetic regions of weak effect. It also increased the phenotypic variance explained for most of the phenotypes considered. It revealed dynamic QTLs with transitory, increasing or decreasing effects over time. By considering both the temporal and genetic multivariate structures in a single statistical model, we increased our understanding of the genetic architecture of complex traits notably by reducing the issue of missing heritability. More broadly, this work raises the foundation for extended applications in functional genomics, evolutionary ecology, and crop breeding programs, in which time-related phenotypic plasticity remains crucial for predicting and selecting key quantitative complex traits. Key messageBy capturing the genetic factors influencing the time related phenotypic plasticity, our approach contributes to a deeper understanding of the dynamic nature of genotype-phenotype relationships.

The actin polymerization machinery, comprising the ARP2/3 complex and its activators, the WASP family proteins, has been implicated in regulating a broad spectrum of nuclear processes, such as transcriptional regulation and nuclear organization. Here, using clustered protocadherin (cPcdh) and {beta}-globin genes as model systems, we showed that WAVE2, a member of the WASP family, regulates chromatin organization by maintaining heterochromatin dynamics. Specifically, by CRISPR DNA-fragment editing, in conjunction with integrated analyses of ChIP-seq, MeDIP-seq, ATAC-seq, 4C-seq, and RNA-seq, we showed that deposition of H3K9me3, a key heterochromatin mark, is significantly decreased at the cPcdh locus upon WAVE2 deletion, concurrent with aberrant accumulation of CTCF/cohesin complex at promoter regions and spatial reorganization of chromatin architecture around nucleolus. In addition, REST/NRSF exerts a similar heterochromatindependent effect on the cPcdh locus. Finally, genetic and genomic data showed that WAVE2 regulates {beta}-globin gene expression by maintaining heterochromatin status. Together our data suggested that WAVE2 and REST/NRSF regulate clustered gene expression in a heterochromatin-dependent manner.

Text abstractsLipid homeostasis is essential for organismal physiology, and its disruption contributes to metabolic disorders. Using an unbiased genetic modifier screen in Drosophila, we identified GAR1, a core component of the box H/ACA small nucleolar ribonucleoprotein complex, as a pivotal regulator of systemic lipid storage. We show that the H/ACA snoRNP complex is essential for maintaining lipid droplet morphology in adipose tissue and preventing ectopic fat accumulation. Moreover, null mutants of Gar1 or Dkc1 exhibit severe developmental defects, including reduced body size and larval lethality. RNA-seq analysis revealed that Gar1 dysfunction triggered widespread alternative splicing defects, specifically targeting key transcripts within the insulin signaling cascade, including chico, Pi3K92E, sgg, and Lip4. Furthermore, knockdown of Gar1 impaired insulin signaling, as evidenced by the reduced membrane localization of the tGPH fluorescence. Genetic epistasis further positions GAR1 upstream of the lin-28/foxo axis, as knocking down lin-28 or foxo fully rescues the lipometabolic defects in GAR1-deficient animals. These findings reveal a previously unrecognized link between the snoRNP machinery and metabolic process, establishing the box H/ACA complex as an important coordinator that integrates RNA processing with insulin-mediated nutrient sensing to ensure developmental and lipid homeostasis. Article summaryLipid metabolism is tightly controlled by multiple factors. To find new regulators, the authors performed a genetic screen and identified a small nucleolar protein GAR1 participate in fat storage and larval development. They demonstrated a critical role of box H/ACA snoRNP complex in modulating alternative splicing and balancing insulin cascade. Blocking two insulin-related genes reversed the lipid defects caused by Gar1 loss. These findings revealed the box H/ACA complex integrates RNA processing with insulin-mediated nutrient sensing to ensure developmental and lipid homeostasis, offering a perspective for understanding the metabolic regulation network.

Gene amplification is a potent driver of evolution and is thought to contribute to genetic diseases, including cancer. The yeast Saccharomyces cerevisiae is a powerful organism for understanding amplification mechanisms. When yeast is grown long term in sulfate-limiting chemostats, amplification of the gene that encodes the primary sulfate transporter, SUL1, is a common outcome. Here we describe a form of SUL1 amplification in which multiple copies of the right terminal region of chromosome II are appended in tandem to a native telomere. We find this form of amplicon when we delete the origin of replication next to SUL1 or delete a variety of genes involved in DNA metabolism. It is the only form of amplification found in a yku70{Delta} mutant suggesting that unprotected telomeres are involved. We propose that these terminal addition events occur when the unprotected 3 G1-3T telomeric sequence invades a short ([~]7 bp) internal telomere sequence (ITS) to begin a form of microhomology-mediated break-induced replication (mmBIR) that has been documented in type-I survivors of telomerase mutants. In addition to amplification of the right end of chromosome II we also find that telomeres containing the sub-telomeric repeat Y experience similar tandem amplification events and show that their formation is reduced in a pol32{Delta} mutant, a gene required for mmBIR. Within individual amplicons the ITSs and Ys are nearly identical, suggesting that the multiple copies of the amplified region are generated in a single mmBIR event that we describe as pseudo-rolling circle mmBIR. A similar amplification event at the P-telomere of human chromosome 18 has four copies of a [~]54 kb region separated by ITSs of nearly identical size. This finding suggests that these additional copies of the terminal fragment of human chromosome 18 arose by the same pseudo-rolling circle mechanism, perhaps during a period of telomeric stress. AUTHOR SUMMARYThe human genome is peppered with duplicates (or higher numbers) of segments that are located at sites both nearby and distant from the original, ancestral segments. These Copy Number Variants, or CNVs, appear to be highly variable among different individuals and are being examined with great interest as potential loci associated with genetic disease. Experimentally determining how these CNVs arise and become distributed across the genome is nearly impossible using humans. We are using budding yeast as the model organism to explore mechanisms of gene amplification. In this work we show that by destabilizing the ends of yeast chromosomes (telomeres) or by interfering with genes involved in the replication, repair, or recombination of DNA results in a specific form of segmental copy number increase that is initiated at telomeres. We propose that a telomere invades an internal chromosome site and sets up a pseudo-circular template for conservative DNA replication. The outcome is a chromosome with multiple, identical copies of a chromosome end arranged in tandem. We believe that it is also a major mechanism used by cells to repair telomeres that have become eroded during aging.

Eukaryotes typically maintain telomere length within a defined range. While short telomeres are known to activate DNA damage responses and limit cell proliferation, long telomeres are associated with extended proliferative capacity. The broader cellular consequences of long telomeres are comparatively less well understood. In budding yeast Saccharomyces cerevisiae, long telomeres have been shown to influence gene expression at specific loci, but whether long telomeres affect transcription genome-wide has not been reported. Here, we analysed transcriptomes in a lineage that inherited long telomeres (originally due to a rif2{Delta} mutation). Transcriptomes were assessed over two rounds of mitosis and meiosis in the absence of the rif2{Delta} mutation. We show that strains with long telomeres exhibit a distinct gene expression profile, including upregulation of membrane transporters and downregulation of a smaller subset of genes. Both up- and down-regulated genes were distributed across the genome, arguing against a purely telomere-proximal effect on gene expression. Affected genes were enriched for Rap1 binding sites, consistent with a model in which long telomeres sequester telomere-associated transcriptional regulators, such as Rap1, and thereby affect gene expression at non-telomeric binding sites for these regulators. Accordingly, the magnitude of transcriptional changes was greatest in strains with the longest telomeres. Together, our findings demonstrate that long telomeres induce a genome-wide transcriptional response that can accompany inherited long telomeres across generations. Similar effects of long telomeres are likely to occur in other eukaryotes, including humans, where long telomeres are associated with disease. Article summaryTelomeres protect chromosome ends, and their length is tightly regulated. While short telomeres are known to be harmful, the effects of long telomeres are less well understood. Using budding yeast, we show that inherited long telomeres alter the expression of dozens of genes across the genome, particularly membrane transporters. These changes are consistent with a model in which long telomeres sequester regulatory proteins away from other loci. Our findings may have broader implications in more complex organisms, including humans.

Microexons are short alternative exons up to 51 nucleotides long that are highly enriched in neuronal genes. Their dysregulation has been linked to human neurodevelopmental disorders, including autism spectrum disorders. In the nematode Caenorhabditis elegans, global regulation of microexons is also critical for proper development. Here we show that microexon alternative splicing (AS) changes between C. elegans larval and adult stages and that microexon inclusion is differentially regulated among distinct neuronal types. Consistently with previous evidence that C. elegans splicing is regulated in response to environmental stimuli, we found here that specific microexons are modulated upon food availability. Both the inclusion levels of these microexons and the feeding behavior seem to depend on the DNA topology, which may affect transcription dynamics, as revealed by the effects of the topoisomerase I (TOP1) inhibitor, camptothecin (CPT). CPT treatment alters responses related to food availability such as speed reduction and exploration. Furthermore, animals carrying a mutation in the global regulator of microexon splicing prp-40 exhibit altered food preference, independently demonstrating that disruption of microexon AS has important consequences on animal behavior.

Cyanobacteria utilize type IV pili for many behavioural responses, such as phototaxis, aggregation, floating, and DNA uptake. Type IV pilus-dependent functions are regulated by the nucleotide second messengers, c-di-GMP and cAMP. In this study, we investigated the role of a recently identified c-di-GMP receptor (CdgR) in cyanobacteria that harbours a ComFB domain. ComFB-domain proteins are widespread in cyanobacteria and are also present in heterotrophic bacteria. We demonstrated that the CdgR homolog from the cyanobacterium Synechocystis sp. PCC 6803, a model organism for studying type IV pilus-dependent functions, specifically binds to c-di-GMP. Genetic and phenotypic analyses revealed that Synechocystis CdgR is involved in phototactic motility and natural competence. Inactivation of cdgR resulted in altered expression of specific sets of minor pilins, which are essential for motility or natural competence. We identified interactions between CdgR and the CRP-family transcription factors, SyCRP1 and SyCRP2. Disruption of these CdgR-SyCRP1 and CdgR/SyCRP2 complexes is initiated by elevated c-di-GMP levels. Moreover, the assembly and stability of these complexes are influenced by other cyclic nucleotides, such as cAMP and c-di-AMP. These observed interactions imply a complex regulatory mechanism by which CdgR influences gene expression in response to cyclic nucleotide messenger signalling, particularly c-di-GMP. The present findings highlight the importance of CdgR in c-di-GMP signalling and its role in regulating type IV pilus-dependent functions in Synechocystis. The modulation of the expression of specific minor pilin genes by CdgR, through interactions with the transcription factors SyCRP1 and SyCRP2, contributes to the establishment of multiple type IV pilus functions and adaptive behaviours of cyanobacteria.

Adaptive mutations, or mutations that confer a fitness benefit, can leave behind distinct signals in genetic data. Computational methods have improved the localization of adaptive mutations in genetic samples using a range of statistical and machine learning classification techniques. However, these methods miss the opportunity to jointly integrate statistics at both the site and window-based level, thus failing to harness all available statistical evidence to detect selection. Our method, WINDEX, combines these different resolutions of statistics to improve the detection of adaptive mutations among hitchhiking signals. Our model simultaneously integrates emissions at different resolutions by defining site-based and window-based latent states corresponding to neutral, linked, and sweep regions, with the site-based states and transition models nested within the window-based states. Using evolutionary simulations with varying selection parameters, we validate the ability of WINDEX to classify positive selective sweeps. Using data from the 1000 Genomes Project, we show that WINDEX is able to identify regions harboring signals of selective sweeps, and provides improved localization within those regions over existing methods. In addition, using WINDEX genome-wide allows for estimation of the proportion of whole genomes that are under positive selective pressures; our estimates of between 9.7-10.5% across different populations provide support for other preliminary estimates of these quantities. Author summaryPopulation geneticists often seek evidence for positive selective sweeps, or an evolutionary event in which a beneficial allele increases in frequency over time in a population, resulting in increased fitness of the individuals that have said allele. Positive selective sweeps, however, are difficult to detect due to varying patterns of linkage disequilibrium (LD), or the nonrandom association of alleles, and detecting these signals reliably among differing LD structures remains a challenge in the field. In this work, we present WINDEX, a probabilistic framework designed to leverage signals of positive selective sweeps at both the site- and window-levels in the form of a hierarchical hidden Markov model (HHMM), to localize regions of positive selective sweeps in aligned haplotype data. We validate WINDEX in evolutionary simulations over varying positive selective sweep scenarios, showcasing the improved resolution that the HHMM structure provides. We apply WINDEX in comparative genomic scans of canonical sites of positive selection as well as whole-genome scans to demonstrate the tools power in localizing functionally-validated signals of selection and to offer insights into the proportion of the human genome currently under positive selective pressures. WINDEX is publicly available and easy to apply to many cases of human genetic data.

Ehmt2 is a key H3K9 methyltransferase that regulates genome silencing and structural integrity during animal development. In addition to this canonical function, Ehmt2 has also been implicated in neural tissues mediating both direct and indirect transcriptional activation, and exon splicing, to facilitate proper neural cell differentiation and survival. Several germline loss-of-function animal models have been developed showing both conserved and divergent phenotypes that range from embryonic lethality to behavioural deficits in adult, fertile animals. Here, we generated the first maternal-zygotic ehmt2 loss of function mutant in zebrafish using CRISPR-Cas9 mutagenesis. An assessment of the pattern of H3K9 methylation in mutant embryos by ChIP-seq indicates that there are aberrant levels of this repressive mark, including reduction in discrete 5 non-coding regions of genes, but with no significant change in the overall pattern distribution of these marks across the genome. Global transcriptome and morphological analyses demonstrated that mutant embryos displayed greater variation in the timing of developmental progression that is, on average, slower compared to controls. Despite this, mutant embryos ultimately survive and are fertile. Through examination of progenitor cell dynamics and gene expression profiles, we found that the delay in embryonic development was associated with longer rates of S-M phases of the progenitor cell cycle in mutants leading to deficits in tissue growth. Finally, our data suggest a robust network of epigenetic regulators can potentially compensate for Ehmt2 loss of function and permit embryonic development and survival in ehmt2 mutant zebrafish. Our work establishes a zebrafish ehmt2 loss of function model that will facilitate examination of the complex and varied roles of Ehmt2 in vertebrate development.

Proper maintenance of gene expression in response to mutations or environmental fluctuations is critical for cell development and survival. Recently, a novel genetic compensation mechanism was described wherein mutant mRNA decay triggers increased transcription of paralogous genes. This effect was reported for several genes, including {beta}-actin (Actb) in mouse embryonic stem cells, where Actb mRNA with a premature termination codon enhances transcription of its paralog, {gamma}-actin (Actg1), and partially rescues cytoskeletal defects. Here we show that, in both mouse and human embryonic stem cells, mutations in the ACTB gene, regardless of mutant mRNA expression, trigger genetic compensation. Furthermore, transgenic expression of mutant ACTB mRNA with a premature stop codon fails to induce genetic compensation. Depletion of the SRF or MRTF-A transcription factors, which are known to increase ACTB transcription in response to low ACTB protein levels, diminishes the genetic compensation response in ACTB mutants. These results suggest that genetic compensation in ACTB mutants is primarily mediated by a transcriptional feedback loop via SRF/MRTF-A, independent of the expression or degradation of mutant ACTB mRNA.

Ribosomes are increasingly recognized as heterogeneous regulators of gene expression, yet how ribosomal protein paralogs interface with nutrient signaling remains poorly understood. In Saccharomyces cerevisiae, ribosomal protein gene expression is governed by duplicated gene pairs, many of which exhibit functional divergence despite high sequence identity. A central regulator of ribosome biogenesis and translational control is the Target of Rapamycin (TOR) pathway, which integrates nutrient signals to modulate growth, stress adaptation, and lifespan. Target of Rapamycin Complex 1 (TORC1) influences ribosome activity by phosphorylating ribosomal protein S6 (Rps6), a modification that links nutrient availability to translational output. Here, we investigated the functional divergence of the Rpl12 ribosomal stalk protein paralogs RPL12a and RPL12b and found that rpl12b{Delta} produces phenotypes consistent with reduced TOR activity, including decreased Rps6 phosphorylation, G2/M cell-cycle accumulation, and significant extension of chronological lifespan as compared to rpl12a{Delta} and wildtype strain. Multi-omics analyses further indicate translational and metabolic reprogramming consistent with activation of a stress-adaptive program associated with Gcn4. Importantly, loss of RPL12b also reduces levels of the ribosome preservation factor Stm1, a TORC1-regulated protein required for stabilization of 80S ribosomes under stress. This finding links ribosomal stalk composition to ribosome stability and nutrient-responsive signaling. Together, our results demonstrate that ribosomal paralog specialization provides an additional regulatory layer connecting translation, TOR signaling, and cellular longevity.

Yeast responding to acute stress reallocate cellular resources, in part via the Environmental Stress Response (ESR) that induces stress-defense genes while repressing ribosome-biogenesis and growth genes. The purpose and regulation of coordinated induction and repression is incompletely understood, but both responses are influenced by ESR transcription factors Msn2 and Msn4 (Msn2/4). Here we used single-cell microscopy and transcriptomic analysis to investigate the role of upstream regulator Pde2 in ESR regulation and post-stress fitness. Loss of PDE2 weakened and shortened Msn2 activation following salt stress and produced muted induction of Msn2/4 targets, similar to a msn2{triangleup}msn4{triangleup} strain. In contrast, Pde2 had at most a minor impact on ESR repressor Dot6, yet was important for repression of its targets beyond Msn2/4 influence. Consistent with our recent resource-reallocation model, pde2{triangleup} cells had normal or faster post-stress growth rates, despite weaker activation of the ESR. We discuss implications for ESR regulation and function.